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Abstract 

Subthreshold kaon (K° and K+) production in neutron-rich nuclear reactions to probe the 
O 1 kaon-nucleon interaction in nuclear medium and to extract the isospin dependent part of the 
, ^, 1 nuclear equation of state at high-baryon densities, is investigated within an isospin and momentum 

,__i ■ dependent transport model. A repulsive kaon-nucleon potential is implemented in the model 

> 

[^ through reproducing the flow data in heavy-ion collisions, which enhances the energetic kaon 

On; 

OO , emission squeezed out in the reaction zone and leads to a variation of the high- momentum spectrum 

t;^ ' of the K^/K"*" yields. It is found that the stiffness of nuclear symmetry energy plays a significant 

cn ' role on the isospin ratio with decreasing the incident energy and a hard symmetry energy has a 

IJ . larger value of the K^/K"*" ratio in the domain of subthreshold energies. 
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Kaon production in relativistic heavy-ion collisions has been investigated as a useful tool to 
constrain the high-density information of isospin symmetric nuclear equation of state (EoS) both 
experimentally [H [2] and theoretically [31 HI O El [7] . Kaons (K^ and K^) as a probe of EoS are 
produced in the high-density domain without subsequent reabsorption in nuclear medium. The 
available experimental data already favored a soft EoS at high baryon densities. The K^ / K^ ratio 
was proposed as a sensitive probe to extract the high-density behavior of the nuclear symmetry 
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energy (isospin asymmetric part of EoS) [HIS], which is poorly known up to now but has an impor- 
tant apphcation in astrophysics, such as the structure of neutron star, the coohng of protoneutron 
stars, the nucleosynthesis during supernova explosion of massive stars etc [10]. Produced kaons in 
heavy-ion collisions can be easily deviated by surrounding nucleons in the dynamical evolutions. 
Consequently, the spectrum of the K^ / K^ yields is to be modified and the constraint of the 
density dependence of symmetry energy from heavy-ion collisions is also influenced. 

In this work, kaon dynamics in heavy-ion collisions is investigated with an isospin and mo- 
mentum dependent transport model (Lanzhou quantum molecular dynamics (LQMD)), in which 
strangeness production is contributed from channels of baryon-baryon and pion-baryon collisions 
[7] . The model has been developed to treat the nuclear dynamics at near Coulomb barrier energies 
and also to describe the capture of two heavy colliding nuclides to form a superheavy nucleus [TTj. 
Further improvements of the LQMD model have been performed in order to investigate the dy- 
namics of pion and strangeness productions in heavy-ion collisions and to extract the information 
of isospin asymmetric EoS at supra-saturation densities [3, [12]. The momentum dependence of 
the symmetry potential was also implemented in the model, which results in an isospin splitting 
of proton and neutron effective mass in nuclear medium [13]. We have included the resonances 
(A(1232), N*(1440), N*(1535)), hyperons (A, S) and mesons (vr, K, t]) in hadron-hadron collisions 
and the decays of resonances for treating heavy-ion collisions in the region of lA GeV energies. 

In the LQMD model, the time evolutions of the baryons (nucleons and resonances), hyperons 
and mesons in reaction system under a self-consistently generated mean-field are governed by 
Hamilton's equations of motion, which read as 

OH . dH 

The Hamiltonian of baryons consists of the relativistic energy, the effective interaction potential 
and the momentum dependent part as follows: 
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Here the Pj and rrii represent the momentum and the mass of the baryons. The effective interaction 
potential f/j„t is composed of the Coulomb interaction and the local potential [13]. The local 
interaction potential is derived from the Skyrme energy-density functional as the form of Uioc = 
J Vioc{p{'r))dr. The energy- density functional reads 
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where the p„, Pp and p = Pn + Pp are the neutron, proton and total densities, respectively, and the 
6 = {pn — Pp)/{pn + Pp) being the isospin asymmetry. The coefficients a, (3, 7, gr, gsur, glur and 
po are set to be the values of -296.6 MeV, 197 MeV, 1.143, 9.9 MeV, 22.9 MeV fm^, -2.7 MeV fm^ 
and 0.16 fm~'^, respectively. A compression modulus of K=230 MeV for isospin symmetric nuclear 
matter is produced with these parameters. A Skyrme-type momentum-dependent potential is used 
in the LQMD model [131 

Umom= Y~ X! ^Cr,r'^T,TA',rj j j dpdp' drfi{r, p,t) 

x[ln(e(p-p')^ + l)f/,(r,p',t). (4) 

Here Cr^r = Cmomi^ + x), Cry = Cmorni^ — x) {j ^ t') and the isospin symbols T(r') represent 
proton or neutron. The parameters Cmom and e was determined by fitting the real part of optical 
potential as a function of incident energy from the proton-nucleus elastic scattering data. In the 
calculation, we take the values of 1.76 MeV, 500 c^/GeV^ for the Cmom and e, respectively, which 
result in the effective mass m*/?7i=0.75 in nuclear medium at saturation density for symmetric 
nuclear matter. The parameter x as the strength of the isospin splitting with the value of -0.65 
is taken in this work, which has the mass splitting of m* > vn* in nuclear medium. The effect of 
the momentum dependence of the symmetry potential in heavy-ion collisions was also investigated 
with the isospin-dependent Boltzmann Uehling Uhlenbeck transport model. The same conclusions 
are found [H]. 

The symmetry energy is composed of three parts, namely the kinetic energy from fermionic 
motion, the local density-dependent interaction and the momentum-dependent potential as 
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The local part is adjusted to mimic predictions of the symmetry energy calculated by microscopical 
or phenomenological many-body theories and has two-type forms as follows: 

e':Up) = lc,ymip/poy% (6) 



and 



Es°ym{p) = asym{p/po) + bsyni{p/po) ■ (7) 



The parameters Cgym, cisym and hsym are taken as the values of 52.5 MeV, 43 MeV, -16.75 MeV. 
The values of 7^=0.5, 1., 2. have the soft, linear and hard symmetry energy with baryon density, 
respectively, and the Eq. (7) gives a supersoft symmetry energy, which cover the largely uncertain 



of nuclear symmetry energy, particularly at supra-saturation densities. All cases cross at saturation 
density with the value of 31.5 MeV. We chose two typical variations with baryon density, i.e., hard 
and supersoft symmetry energies in the domain of high densities. 

The hyperon mean-field potential is constructed on the basis of the light-quark counting rule. 
The self energies of hyperons are assumed to be two thirds of that experienced by nucleons. Thus, 
the in-medium dispersion relation reads 



a;(p„ pO = ^{mn + ^If + p? + S^ (8) 

with Ef^ = 2E;5/3 and Ey = 2'Ly /?>, which leads to the optical potential at the saturation 
density being the value of -32 MeV. The evolution of mesons (here mainly pions and kaons) is also 
determined by the Hamiltonian, which is given by 

Hm =E(^.^°^^ + ^(p.a)). (9) 

Here the Coulomb interaction is given by 

^Coul ^ g ^^ (^Q) 

where the Nm and A''^ are the total numbers of mesons and baryons including charged resonances. 
We consider two scenarios for kaon (antikaon) propagation in nuclear medium, one with and one 
without medium modification. From the chiral Lagrangian the kaon and antikaon energy in the 
nuclear medium can be written as [U [T5] 

^KijPi,Pi) 
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m-^ + Vi- a-Kpi + ipKpi) \ - bxpi, (12) 

respectively. Here the Bk = 3/(8/^) ^0.32 GeVfm^, the ax and aj^ are 0.18 GeV^fm^ and 0.3 
GeV^fm^, respectively, which result in the strengths of repulsive kaon-nucleon (KN) potential and 
of attractive antikaon-nucleon potential with the values of 25.5 MeV and -96.8 MeV at saturation 
baryon density, respectively. The values of m*j^/mK=^-05 and 171^/171-1^=0.8 at normal baryon 
density are concluded with the parameters. The effective mass is used to calculate the threshold 
energy for kaon and antikaon production, e.g., the threshold energy in the pion-baryon collisions 

The scattering in two-particle collisions is performed by using a Monte Carlo procedure, in 
which the probability to be a channel in a collision is calculated by its contribution of the channel 



cross section to the total cross section. The primary products in nucleon-nucleon (NN) colhsions 
in the region of lA GeV energies are the resonances of A (1232), A^*(1440), A^*(1535) and the 
pions. We have included the reaction channels as follows: 

NN ^ NA, NN ^ NN*, NN ^ AA, 

A ^ A^vr, A^* ^ Ntt, NN ^ A^A^7r(s - state), 

iV*(1535) -^ Nr]. (13) 

At the considered energies, there are mostly A resonances which disintegrate into a tt and a 
nucleon in the evolutions. However, the A^* yet gives considerable contribution to the energetic 
pion yields. The energy and momentum-dependent decay widths are used in the model for the 
resonances of A(1232) and A^*(1440) [12]. We have taken a constant width of 1=150 MeV for the 
A^*(1535) decay. 

The strangeness is created as the secondary products in inelastic hadron-hadron collisions. We 
included the channels as follows: 

BB -^ BYK, BB -^ BBKK, Bn -^ YK, Bn -^ NKK, 

Ytx -^ BH, B'K -^ Ytx, YN -^ 'KNN. (14) 

Here the B stands for (N, A, N*) and Y(A, S), K(K°, K+) and 'K{W , K"). The elastic scattering 
between strangeness and baryons are considered through the channels of KB — )■ KB, YB — )■ YB 
and KB — t- KB. The charge-exchange reactions between the KN — )■ KN and YN — )■ YN 
channels are included by using the same cross sections with the elastic scattering, such as K^-p — )■ 
K^n, K^n — )■ K^p etc. Correction of effective mass of kaons in nuclear medium on the elementary 
cross section is considered through the threshold energy, which results in the reduction of kaon and 
the enhancement of anti-kaon yields in heavy-ion collisions. Shown in Fig. 1 is a comparison of the 
total kaon yields produced in the central ^^^Au+^^^Au collisions with and without KN potential 
for the hard (left window) and supersoft (supersoft window) symmetry energies, respectively. The 
multiplicity of K" is multiplied by 10 times. Inclusion of the KN potential in the model leads to 
about 30% reduction of the total kaon yields in the subthreshold domain. Once kaons are produced 
in the compression stage, the subsequent reabsorption by the surrounding baryons almost does not 
take place although the distributions of kaons in phase space can be deviated by the in-medium 
potential. 

To investigate kaon production in momentum space and its correlation to the collision geometry 
and to the KN potential, we computed the rapidity and transverse momentum distributions of 



K° and K"*" in the near central ^^^Au+^^^Au reaction (b=l fm) with a hard symmetry energy as 
shown in Fig. 2 in the left windows. The spectra of inclusive invariant cross sections in the mid- 
rapidity region for the K"*" production as measured by the KaoS collaboration [1] are compared 
with the LQMD calculations in the right windows. A nice agreement between the experimental 
data and the calculations with inclusion of the KN potential can be seen. The difference of K*^ 
and K"*" and the influence of nuclear symmetry energy on the kaon inclusive spectra are very 
small. One notices that the repulsive KN potential reduces the kaon yields in the mid-rapidity 
region. A broad rapidity distribution and a flat transverse momentum spectrum are found for the 
case with inclusion of the KN potential. It is caused from the fact that the repulsive potential 
enhances the energetic kaon emission and reduces the kaon yields owing to the threshold effect. 
The isospin effects are pronounced in the domain of mid-rapidity, in particular for the case without 
the KN potential, where kaons are squeezed out from the reaction zone and produced mainly at 
supra-saturation densities formed during the compression stage of two colliding nuclides. The 
azimuthal anisotropy of kaon emission is presented from the transverse flow as shown in Fig. 3. 
The same trends with proton flows but with smaller values are observed for the case without the 
KN potential. However, the repulsive potential even leads to the appearance of anti-fiows. It 
should be noticed that the nuclear Lorentz force is missed in this work, which plays a significant 
role on the kaon flows in heavy- ion collisions [TB]. 

Kaons are produced at the early stage in heavy-ion collision and promptly emitted after produc- 
tion, which can get directly the information of high-density phase diagram [7] . More pronounced 
effects of the KN potential and the stiffness of symmetry energy on kaon production can be ob- 
served from the spectrum of isospin ratio. Shown in Fig. 4 is the influence of the KN potential on 
the transverse momentum distribution of the K^/K"*" ratio for the hard symmetry energy in the 
left panel and the supersoft case in the right window. In the both cases the KN potential plays 
an important role at high transverse momenta and a flat spectrum appears after inclusion of the 
potential. The results will be helpful for experimental measurements to extract the in-medium 
kaon potential. The isospin effects appear at deep subthreshold energies as shown in Fig. 5. 
The in-medium potential slightly changes the K^/K"*" value because of its influence on the kaon 
propagation and also on the charge-exchange reactions. At the considered energies, the channel of 
A^A — i- NYK contributes the main part for the kaon yields due to the larger production cross sec- 
tions and the higher invariant energy, and the A^A^ — )■ NYK as well as iiN — )■ YK have about one 
third contributions. One notices that a hard symmetry energy always has the larger values of the 
isospin ratios than the supersoft case in the domain of subthreshold energies {Eth{K)=l.'b'8: GeV). 



This is caused from the enhanced production of A^ resonances, i.e., nn — > pA^, pA^ — t- riKK^ . 
The threshold effect enlarges the isospin effect of the K°/K+ yields [8]. 

A comparison to experimental data is performed from the double ratio of K^/K"'' taken from 
two systems of ^^Ru+^^Ru and ^^Zr+^^Zr as shown in Fig. 6. It should be noted that the 
influence of symmetry energy on the excitation function of double ratio is very weak over the 
whole energy range. The less sensitivity to the stiffness of symmetry energy is consistent with the 
calculations using the RBUU transport model [IT] . It is well known that the kaon yields increase 
with the number of participate nucleons of colliding nuclides. A larger high-density zone can 
be formed in heavier colliding systems, where kaons are produced in baryon-baryon and meson- 
baryon collisions, and the influence of the stiffness of symmetry energy on the isospin kaons is 
pronounced. Therefore, the isospin ratio in heavy neutron-rich reactions can be a nice probe to 
extract the density dependence of nuclear symmetry energy. 

In summary, kaon dynamics in heavy-ion collisions at near threshold energies has been investi- 
gated by using an isospin- and momentum-dependent transport model (LQMD). It is found that 
the KN potential plays an important role on kaon emission in phase space, in particular, reducing 
the kaon yields in the mid-rapidity region and also at low transverse momenta, but enhancing 
the production at high transverse momenta. The isospin ratio of K^ / K'^ depends on the kaon 
potential at high transverse momenta and a flat spectrum appears after inclusion of the KN po- 
tential. The K^ / K^ ratio of neutron-rich heavy system in the domain of subthreshold energies is 
sensitive to the stiffness of nuclear symmetry energy, which is a promising probe to extract the 
high-density information of symmetry energy through comparison to experimental data. Precise 
measurements on subthreshold kaon production from neutron-rich nuclear collisions are still very 
necessary. 
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Figure 1: Total multiplicities of K'^ and K"*" as a function of incident energy in central ^^^Au+^^^Au 
collisions at the cases of hard (left panel) and supersoft (right panel) symmetry energies. 
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Figure 2: Longitudinal rapidity and transverse momentum distributions of isospin kaons (left 
panels) and the kinetic energy spectra of inclusive invariant cross sections at mid-rapidity (right 
panels) in the ^^"^Au+^^^Au reaction at the incident energy of 1.5^ GeV. The mid-rapidity condi- 
tion is a selection of 6c.m. = 90° ± 10° both for the data and the calculations. 
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Figure 3: Transverse flows of kaons produced in tlie ^^"^Au+^^^Au reaction at tlie incident energy 
of 1.574 GeV with and without the KN potentials as a function of the longitudinal rapidity. Proton 
flows are also implemented for a comparison. 
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Figure 4: Comparison of transverse momentum distributions of the K'^/K"'" ratio with and without 
the KN potentials in central ^^^Au+^^^Au collisions at the incident energy of 1.5^1 GeV. 
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Figure 5: (Color online) Comparison of excitation functions of the K^ / K^ yields for central 
^^^Au+^^^Au collisions for the cases of hard and supersoft symmetry energies. 
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Figure 6: Excitation functions of the double isospin kaon ratios taken from the reactions of 
^^Ru+^^Ru and ^^Zr+^^Zr in the cases of hard and supersoft symmetry energies and compared 
with the experimental data [T7] . 



I.ti 



